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ABSTRACT 


Composite materials are seeing inereased use in struetural applieations beeause of 
their various benefits. When eomposite structures are employed in a water environment, 
their dynamic responses are greatly affected by the fluid m edium. Water density is 
comparable to many composite materials and the effects of fluid-structure interaction on 
dynamic behaviors of composite structures are significant. The effect s of fluid-structure 
interaction include changes of frequency, magnitude, en ergy dissipation, etc., of 
structural characteristics. Hence, it is critical to understand the fluid-structure interaction 
of composite structures subjected to dynamic loading in water environments. This work 
focuses on finding param eters affecting the transient dynam ic responses of com posite 
structures. Coupled fiuid-st ructure interaction analyses of com posite plates are 
conducted numerically, using finite element models, including various parametric studies. 
The results are compared to those of dry structures to identify the role of each parameter. 
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I. INTRODUCTION 


A. BACKGROUND 

Composites are seeing inereas ed use in m aritime, aerospaee and a utomotive 
structures used in both civil and m ilitary applications. Early uses of c omposites were 
limited to secondary structures; however, as knowledge and understanding of mechanical 
characteristics of composites has grown, more primary load-bearing structures have been 
fabricated. In recent years, la rge composite structures have been incorporated into naval 
vessels to increase operational perf ormance wh ile lowering ownership costs [1]. For 
example, carbon-fiber com posite material pr ovides high strength and stiffness while 
maintaining low weight, which in tur n translates to inereas ed fuel econo my or increased 
payload. A further advantage of com poshes over m etals is lower m aintenance and 
resistance to corrosion, making composites very desirable for maritime applications. The 
use of com posites in engineering components has initiated num erous studies to analyze 
structural components fabricated from various com posites rather than traditional m etals. 
While composites provide advantages over m etals, they als o come with com plex and 
challenging engineering problems for analysts and designers [2]. Because the stru ctural 
behavior is im pacted by Flui d S tructure In teraction (FSI), th is work focuses on the 
implications of utiliz ing com posite structu res in m aritime applications below the 
waterline. 

B, LITERATURE SURVEY 

It is critical to assess the structu ral be havior of com posite structures used in 
marine applications beneath the wa terline wher e FSI p lays an im portant ro le on the 
dynamic response and failure of the submerged composite structure. Because composite 
structures are m uch lighter than m etallic structures, the effect of FSI is m uch greater. 
Many polymer composite materials are on ly a few times heavier than water; the refore, 
the added mass effect of the fluid becomes critical. 

Numerous studies hav e exam ined the e ffect of FSI f or m etallic struc tures, 
especially f or underwater explosive loading [3]-[I0]. Som e works are experim ental 
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studies, while others are num erieal work. A few studies exam ined FSI for eom posite 
structures subjected to underwater explosion [11]-[17]. 


C. OBJECTIVES 

This work investigates the effects of the surrounding fluid on dynam ic responses 
of composite structures subjected to a mechanical loading via applied concentrated force, 
uniform pressure and i mpact. The research exam ines several param eters affecting 
transient dynam ic responses of subm erged composite s tructures to iden tify m ajor 
controlling param eters of FSI. Thi s resear eh focuses on com putational m odeling of 
coupled fluid-structure interaction analys es of com posite structures—specifically 
plates—under water for various param etric studies. Results are norm alized to those of 
completely dry structures to illustrate the role of each parameter on FSI. 
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II. COMPUTATIONAL MODEL 


A, MATERI AL SPECIFICATIONS 

The composite material used in this study is an e-glass woven fabric with a plain 
weave fiber architecture and vinyl-ester resin. The com posite has elastic m odulus 17 
GPa, Poisson’s ratio 0.3, and density of 2020 kg/m^. To make a fair comparison between 
dry and wet structures, any potential change of composite material properties associated 
with moisture absorption from water is not considered. The steel used for im pact study 
has elastic modulus 200 GPa, Poi sson’s ratio 0.3, and density 8000 kg/m . For a dry 
structure, i.e., in a ir, there is no spe cific modeling of the air m edium. For m odds that 
examine FSI, the water is modeled with a density of 1000 kg/m , and bulk modulus of 2.2 
GPa, while water viscosity is neglected. 

B, FINITE ELEMENT MODEL DEVELOPEMENT 

As an initial step in studying Fluid Structure Interaction effects, only linear elastic 
behavior is considered in th is study. Solid m aterials are modeled using the Lagrangian- 
based finite element method, while fluid is solved using the Eulerian-based finite element 
method [18]. The com posite plate used in the study is thin (0.002 m thickness), having 
an aspect ratio of at least 150 (length to thickness), and nece ssitates modeling through 
shell elements. 

Due to the thin com posite plate re quiring to b e m odeled with shell e lements, 
coupling between the fiuid/com posite interfaces presented a challenge, as th e interface 
between them needs to be uniquely defined by a volum e or solid clem ents. To have a 
uniquely-defined volume, a stiffened composite plate is used to create a unique volume of 
composite. The stiffened com posite shell structure is composed of top and bottom skin 
plates, coupled through vertical stiffeners. Each skin plate is 0.3 m x 0.3 m and 0.002 m 
thick and is m odeled with e-glass com posite. The stiffeners are m odeled of the sam e 
composite material with the same thickness of the skin plates. Their sizes are 0.3 m long 
and 0.01 m tall, and spaced every 0.05 m apart. The spacing between nodes of composite 
model was 1 cm , such that the 0.3 m skin pi ates have a 30 by 30 m esh. The stiffened 

3 



composite plate is depieted in Figure 1, with the lines denoting the loeations of stiffeners. 
With this stiffened eomposite plate, FSI ean be investigated by eomparing three different 
eases: 1) completely dry, 2) two-sides wet, and 3) one-side wet. 



Figure 1. Stiffened Composite Plate Strueture 

Various parametrie studies, including boundary conditions and loading types, are 
examined to investigate FSI effeets. The edges of the stiffened eomposite plate use either 
a clamped or sim ply-supported boundary eondi tion. The plate is s ubjected to eonstant 
applied force at the een ter of the top skin plate, equivalent pressure loading over the 
surface of the top plate, or impaet loading at center of the top plate from a steel projeetile 
at various initial veloeities. 

The Finite Elem ent Models (FEM ) were eonstructed in PATRAN and solved 
numerieally using DYTRAN. The eom putations were run using a HPC cluster system . 
The computational time required to perform 0.05 second transient solutions varied from 
approximately 5 minutes for the dry strueture to as mueh as 40 hours for the one-side wet 
structural model. The dry ease struetural model has 2,220 elem ents and the wet m odels 
have up to 30,000 elem ents. The geometry used to define the eomposite material uses a 
Eagrangian-based quadrilateral shape for defining the shell elements. The geometry used 
to define the Eulerian-based fluid is eomposed of hexagonal solid elements. 
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1, Dry Structure 

The reference case throughout the study is a completely dry structure using only 
the composite plate constructed as describe d previously. No sp ecific m odeling of air 
surrounding, and within , the void spaces of th e stiffened com posite plate structure is 
accounted for, due to its negligible effects. The dry structure dynam ic response is used 
for normalization with other cases to show the effects of FSI. 

2. Two-Sides Wet Structure 

A two-sides wet structure is used to exam ine the influence of fluid (water) on the 
response. It is m odeled with the stiffened composite plat e embedded within a cube of 
water. The surrounding fluid dom ain is much greater than the composite plate structure 
with a two to one ratio of largest dimension. Addition ally, the non-reflectiv e boundary 
condition is applied to the outside fluid boundary. Although there may be some reflected 
waves from the non-reflective boun dary due to im perfect boundary condition, th e time 
period of interest for structural response is too short to include the effe cts of reflected 
waves. 


3, One-Side Wet Structure 

A one-side wet structure is used to simulate a condition in which fluid is on one 
side of the plate while air is on the other, such as would be encountered in construction of 
a ship hull with com posite p lates. To create an air sp ace on one side of the stiffened 
plate, five additional rigid composite sides are added below the stiffened plate. The sides 
are rigid sh ells composed of the sam e composite m aterial and form the volum e to be 
coupled with the surrounding fluid. The one-sid e wet structure is depicted in Figure 2. 
The air volume between the bottom of the stiffened plate and the bottom of the rigid box 
is 0.01m in height. 
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(a) Box made of a stiffened composite plate and five rigid sides 



X 

z 

(b) Composite box inside water 

Figure 2. Stiffened Composite Plate Supported by Rigid Box used in One-Side Wet 

Case 
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III. PARAMETRIC STUDIES USING COMPUTATIONAL MODEL 


A, TYPE OF BOUNDARY CONDITION 

Two different boundary eonditions are applie d to the stiffened eom posite plate, 
elamped or sim ply supported. In reality it is diffi eult at best to aehieve a perfeetly 
elamped boundary eond ition, and th us aetual boundaries are a m ixture of elarn ped and 
simply supported. To bound the dynam ic response of eom posite plate, both boundaries 
are app lied individu ally to determ ine a ny differenee between FSI effeets. Any 
experimental work done in conjunction with this study will have im perfectly clamped 
boundaries, and thus the behavior will be a mixture of both boundary co nditions. These 
numerical models can be used to understand the differences. 

B, APPLIED LOADING TYPE 

The basis for this study uses an applied concentrated force of lOOON at the center 
of the top skin plate to observe the dynam ic response and determ ine the FSI. 

Additionally, an equivalent pressure to the concentrated force is also examined to reveal 
any differences in response from loading methods. Finally, im pulse type loads are 
imparted to the com posite plate using steel projectiles. The steel projectiles are 0.3 m 
long, and have either a circular or square impact face with area of 1.6129e' m (0.25 in ). 
The steel projectiles start 2 mm above the top skin plate, and are given an initial velocity 
of 1 m/s, 5 m/s or 10 m/s. 

C. PLATE SIZE 

The basic stiffened composite plate used in this num erical study consists of a 
0.3 m by 0.3 m skin plate. A larger 0.5 m by 0.5 m skin plate model is also examined, so 
the differences in FSI can be examined from increased spacing between supports. 

D. PLATE SHAPE 

The basis for this study is the standard 0.3 m by 0.3 m square stiffened plate. To 

examine the impact of plate shape on the dynamic response and FSI, an equivalent area 

rectangular shaped plate is also modeled with dimensions of 0.2 m by 0.45 m. 
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E, COMPOSITE MATERIAL PROPERTIES 

Parametric studies are conducted using the basic 0.3 m by 0.3 m stiffened plate to 
examine the effect of composite m aterial properties on FSI and dynam ic response. T he 
composite material is modeled with a nominal density of 2020 kg/m and nominal elastic 
modulus of 1.7e^° Pa. Two different densities, approximately a 50% reduction and 100% 
increase from the nom inal, are used to inve stigate the change in resp onse; specifically, 
the composite densities are 1020 kg/m and 4020 kg/m . Two different elastic m oduli, 
approximately a 50% reduction and 50% in crease from the nom inal, are used to 
investigate the change in response; specili cally, the composite elastic m oduli are 0.7e 
Pa and 2.7e'° Pa. 
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IV. NUMERICAL STUDY RESULTS AND DISCUSSION 


A, METHODOLOGY 

The basic methodology used to determine the difference in dynamic behavior of a 
stiffened composite pla te is to normaliz e the tw o-sides and one-side we t cases with the 
completely dry case. I n this m anner, the dry case is the base re sponse and tends to 
differentiate the particular changes due to th e FSI. In general, th e base case us ed for 
normalization is the completely dry plate with composite properties of 2020 kg/ m3 for 
density, 1.7e Pa for elastic m odulus with clampe d edges, and a 1000 N concentrated 
force applied at the center of the top skin p late. W hen strains are exam ined the 
normalization is accomplished with respect to the normal x-axis strain. 

This method of norm alization shows the tr ansient variation of various response 
variables; such as displacement of the central node of top skin plate, strain energy and/or 
kinetic energy of the stiffened composite plat e, and stress or stra in at one of three 
locations on the bottom skin plate. The numerical solutions from DYTRAN using shell 
elements only perm it stress to be determ ined. Strains are calculated us ing the s tandard 
stress/strain transf ormation equatio ns. In the com putational m odd, stress (and hence 
strain through transformation equations) is calculated at th e element in center of plate of 
one quadrant (this location is termed ‘center’), at an element half way between the center 
and edge of one quadrant (this location is te rmed ‘side’), and at an elem ent half-way 
between the center and the corner along a di agonal of one quadran t (this location is 
termed ‘quarter’). An exam pie of this sch erne of specific elem ents used to calcu late 
stress/strain is shown in Figure 3 for a 10 by 10 element mesh, although actual composite 
plate mesh is finer. 

The normalized transient responses of di splacement and strain energy typically 
show the same shape and frequency, with onl y minor differences in relative am plitudes, 
and thus can be used interchangeably to demonstrate the behavior of the composite plate. 
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Figure 3. Sample of Element Loeations used to Caleulate Stress/Strain 

B, DYNAMIC RESPONSE OF COMPOSITE PLATE SUBJECTED TO 
CONCENTRATED FORCE AND CLAMPED BOUNDARY 

The baseline stiffened eomposite plate of density 2020 kg/m3 and elastie modulus 
of 1.7e P a with clamped edges and centrally applied concen trated force w ill be 
discussed first. Folio w on sectio ns will exam ine variations in boundary condition , 
loading, size, shape and im pact. Figures 4, 5 and 6 show the response of the 
displacement, strain energy and kinetic energy of the plate respective! y (the dry case is 
used for normalization). 



Figure 4. Normalized Displacement at Center of Top Skin Plate 
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Figure 5. Normalized Strain Energy of Composite Plate 



Figure 6. Normalized Kinetic Energy of Composite Plate 

These figures show the com parison be tween one-side wet and two-sides wet 
structural responses. The ESI with either one-side or tw o-sides wet of the com posite 
structure significantly influences both the m agnitude and frequency of the strain energy 
plot. The oscillating m agnitude and the fir equency are drastically reduced by the ESI 
effect. Two-sides wet ESI results in the lowe st peak energy values and their frequency 
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among the three eases. However, the m agnitude of oseillatory behavior is the least for 
the one-s ide wet stru eture. The figures sh ow the effeets of FSI, with average 
displacement and energy being reduced thr ough the fluid in teraction. Ad ditionally, FSI 
causes a decrease in frequency and m agnitude of structural responses, with significantly 
more rapid damping effects than the dry case. 

The transverse displacem ent plot at the node of the applied force is compared in 
Figure 4 for three different cases. The displace ment response is very sim ilar to that of 
strain energy of Figure 5. The two-sides wet structure has the lowest peak displacem ent 
and frequency, and the one-side w et structur e has the least vibrat ory m otion. It is 
interesting to note that even though the disp lacement characteristics are quite different 
among the three conditions; their respective average values are comparable. 

When average values of the th ree strain energy variations are com pared (Figure 
5), the dry structure has the greatest average value and the two-sides wet structure has the 
smallest value. Furthermore, the two-side s wet structure shows energy dissipation as a 
function of time. 

As the kinetic energy of the stiffened st ructure is compared under three different 
surrounding m edia, as shown in F igure 6, the dry structure show s a very significant 
oscillatory behavior. O n the other hand, the oscillation of kinetic energy is suppressed 
quickly for the wet cases. The kinetic energy of the two-sides wet structure is the lowest. 
The two-sides wet structure displays the fastest decay rate of the kinetic energy. 

The normal and shear strains for th is clamped case, for each of the locations of 
interest (center, side and quarter), are shown in Figure 7, with normalization with respect 
to dry plate .r-axis normal strain. Com parison of the norm al strain along the v-axis also 
indicates reduced strain s for wet structures. Wet structures have very high frequency 
components in the strain response. Howeve r, the base frequencie s of both-side wet 
structures are clearly shown lower than those of the dry structure. Average strain v alues 
are more or less sim ilar even though the dry structure has greater am plitudes of strain 
oscillation. 
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Figure 7. Normal and Shear Strains for Clamped Boundary with Coneentrated Foree 
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C. CLAMP VERSUS SIMPLE SUPPO RT B OUNDARY CONDITION WITH 

CONCENTRATED FORCE 

Comparison of cla mped versus simple boundary shows little difference in 
dynamic response for concentrated force load ing. The displacem ent and strain energy 
plots are shown in Appendix A. The kinetic energy responses shown in Figures 8, 9 and 
10 shows the com parison for the dry, two-sides wet and one-side wet structures 
respectively. There is aim ost no differen ce for t he dr y st ructure; however , t he wet 
structures show slight increase in energy. This is expected due to the increased degree of 
freedom, although the increased energy is not significant. 



Figure 8. Comparison of Kinetic Energy of Dry Structure between Clamped and 

Simple Boundary 



Figure 9. Comparison of Kinetic Energy of Two-sides Wet Structure between 

Clamped and Simple Boundary 
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Figure 10. Comparison of Kinetic Energy of One-side Wet Structure between 

Clamped and Simple Boundary 

Similarly, the strains at the center an d side locations are nearly identical and are 
shown in Appendix A. Of interest are the quarter location strains, which show som e 
variance between the boundary condition types, with the clamp condition having slightly 
higher strains for the dry and wet structures as shown in Figure 11. The increase in strain 
for the clamped boundary was expected due to restricted degree of freedom; however it is 
surprising to be evident at only the quarter location. 

With an applied conce ntrated f orce, th ere is little d ifference betwe en the two 
types of boundary conditions, clamped or simple support. While there is minor increase 
in kinetic energy of the wet cas es for simple support and m inor decrease in strain at the 
quarter location of the com posite plate for the sim pie support, it is not significant. The 
ESI effects are consistent between the two boundary conditions. 
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Figure 11. Normal and Shear Strain Comparison at Quarter Position for Clamped 
versus Simple Boundary with Concentrated Force 
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D. 


CONCENTRATED F ORCE VERSUS P RESSURE LOADING WITH 
CLAMPED BOUNDARY 


Next, the dynamic response of thin composite plate was compared under different 
loading con ditions: con slant concentrated fo rce and equiv alent uniform pressure, each 
with clamped boundary. The basis for com parison is clam ped boundary with constant 
concentrated force of 1000 N applied at center of plate. The equivalent uniform pressure 
loading is determ ined from the concentrated force being uniform ly a pplied over the 
surface of the 0.3m by 0.3m plate, giving a uniform pressure load of 11,111 Pa. 

The comparison for the dry structure under the two loading conditions is shown in 
Appendix B. Under dry conditions, the pre ssure loading versus concentrated force 
increases the amplitude of oscillation for displacement, strain energy, and kinetic energy 
with no shift in frequency. The strain at the center location has increased amplitude but 
lower average strain. T he normal average stra in at the side location is incre ased, while 
the shear strain is com parable between the two loading conditions. The quarter location 
exhibits similar strain behavior for applied force and pressure loading. 

The dry structure is used to norm alize the wet structure responses and the 
displacement response showing the FSI effects are shown in Figure 12. The wet structure 
comparison of strain and kinetic energy for force versus pressure load is shown in Figures 
13 and 14, respectively. 




Figure 12. Wet Structure Displacement Comparison between Force and Pressure 

Loading with Clamped Boundary 
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Figure 13. Wet Structure Strain Energy Comparison between Force and 
Pressure Loading with Clamped Boundary 




Figure 14. Wet Structure Kinetic Energy Comparison between Force and 
Pressure Loading with Clamped Boundary 

Figures 12 to 14 show t he comparison between one-side wet and two-sides wet 

structural responses. The FSI with both one-side and two-sides wet reduces the 

oscillating magnitude and frequency of the response over dry structure. The pressure 

load tends to produce larger am plitude of os cillation than concentrated force, but th e 

average energy is similar, while the mean displacement under pressure load is less. Two- 

sides wet FSI results in the lowest peak energy, peak displacement, and frequency among 

the three cases. However, the m agnitude of oscillatory behavior is the least for the two- 

sides wet structure. The figures sho w the effects of FSI, wi th average displacement and 
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energy being redueed through the fluid intera etion. Additionally, FSI eauses a deerease 
in frequeney and m agnitude of struetural responses with signifieantly m ore rapid 
damping effeets than the dry ease. Of note in Figure 14 is the slower initial response of 
kinetie energy under pressure load. 

The strain behavior for the dry strueture is shown in Appendix B. The pressure 
load viee eoneentrated foree eomparison show inereased amplitude of strain oseillation at 
eenter loeation and redueed average strain as well. The strain response at the side and 
quarter loeations was similar for both normal and shear strains, with exeeption of norm al 
strain at the side loeation ha ving a slightly higher m agnitude under applied pressure than 
applied foree. The strain behavior f or wet struetures, also shown in Appendix B, was 
similar to that of dry, with an inereased amplitude oseillation at eenter loeation, but with 
redueed average strain. The side and quarter loeation strain response of wet struetures 
also follow ed that of dry, with sam e diff erenees of the norm al side loeation s train 
exhibiting higher m agnitude under pressure lo ading. This m eans the eoneentrated foree 
has a greater FSI effeet than the pressure at the side and quarter loeations. 

E. CONCENTRATED F ORCE VERSUS P RESSURE LOADING WITH 

SIMPLE SUPPORT BOUNDARY 

Next the dynam ie response of a thin eom posite plate was eom pared under 
different loading eonditions, eonstant eon eentrated foree and equivalent unif orm 
pressure, with a sim pie support boundary. The basis for com parison is sim pie support 
boundary with constant concentrated force. 

The com parison for the dry structu re under th e two loading conditio ns with 
simple support is shown in Appendix C. Under dry conditions the pressure and 
concentrated force loading have nearly identi cal responses with no di scemable change in 
amplitude or frequency for displacement, strain energy and kinetic energy. The strains at 
the center locations have nearly identical response for applied force and pressure loading. 
The norm al averag e strains at th e side and qu arter lo cation is in creased for press ure 
loading, wh ile the shear strain sho ws higher amplitude of oscillation for the pres sure 
loading condition. 
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The wet structure response com parison with simple boundary for force versus 
pressure load is shown in Figures 15, 16 and 17 for displacem ent, strain energy and 
kinetic energy respectively. 




Figure 15. Wet Structure Displacement Comparison of Force and Pressure Loading 

with Simple Support Boundary 




Figure 16. Wet Structure Strain Energy Comparison of Force and Pressure Loading 

with Simple Support Boundary 
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Figure 17. Wet Structure Kinetic Energy Comparison of Force and Pressure Loading 

with Simple Support Boundary 

The ESI with both one-side and two-sides wet structures reduce the os cillating 
magnitude and frequency of the response over the dry structure. With a simple boundary, 
the pressure and force load track very well with one another, with only minor difference 
in frequency evident in the two-sides wet structure displacem ent and strain energy. 
Unlike the clam ped boundary, there is no dela y in response of kinetic energy with a 
simple boundary for the force and pressure loa d. Again, the two-side s wet ESI results in 
the lowest peak energy, peak displacem ent and frequency among th e three cases. The 
figures sho w the effects of ESI, with aver age displacem ent and energ y being red need 
through the fluid interaction. 

The strain behaviors for the three stru ctures are shown in Appendix C. The 
pressure and force load strains track each other using a sim pie support boundary at the 
center location. The strain respon se at th e side location has sim ilar am plitude of 
oscillation, with the norm al st rains slightly higher under pr essure load. The strain 
behavior at the quarter location is sim ilar for wet structures, although the two-sides wet 
structure has less amplitude, the wet structures overall have approximately equal average 
strain. 

F. SIZE OF COMPOSITE PLATE 

Next the influence of com posite plate size on ESI is exam ined by increas ing the 

size of the square plate from 0.3m to 0.5m on a side. The com parison is m ade using 
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clamped boundary condition with applied concentrated force. The displacement response 
for the three structures is shown in Appendix D, and indicates that increases in plate size 
yield a decrease in frequency, with the tw o-sides wet structure having a substantial 
decrease in frequency. Also, FSI damping is slower as the plate size in creases. Similar 
results are visible in strain and kinetic en ergy response betw een the two sizes of plates 
shown in Figures 18, 19 and 20, for the dry, two-sides wet and one-s ide wet structures 
respectively. 




Figure 18. Dry Structure Strain and Kinetic Energy Comparison for Plate Size 
Variations with Concentrated Force and Clamped Boundary 




Figure 19. Two-sides Wet Structure Strain and Kinetic Energy Comparison for Plate 
Size Variations with Concentrated Force and Clamped Boundary 
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Figure 20. One-side Wet Structure Strain and Kinetic Energy Comparison for Plate 
Size Variations with Concentrated Force and Clamped Boundary 

The larger size plate has lower frequenc y, slower initial response, and slower 
long-term damping. The ESI of the larger pi ate is less and thus has higher am plitude 
oscillations. The dry structure has similar average energies between the two plate sizes, 
while the steady state energy of the larger plate is marginally greater for the one-side wet 
structure. The difference in energy between the two plate sizes is more pronounced in the 
two-sides wet structure, where the kinetic energy clearly shows th e significant delay in 
response due to the damping effect of fluid. 

The comparison of strain between the two sized plates is shown in Appendix D. 
However, there is no clearly identifiable characteristic between the strains with exception 
of some decreased frequency and comparable average normal and shear strains. 

G. SHAPE OF COMPOSITE PLATE 

The influence of composite plate shape on ESI is exam ined next by changing the 

shape of the plate from square to rectangular while m aintaining equivalent area, thus the 

rectangular plate is 0.2 m by 0.45m. As with com parison of plate size, the shape 

comparison is made using clamped boundary conditions with applied concentrated force. 

The displacement response for the three structures is shown in Appendix E, and indicates 

the rectangular shape has increase in frequency and decrease in amplitude of oscillation 

over the square plate of equivalent area, with the average displacem ent of the th ree 

structures (dry, two-sides wet, one-side wet) slightly greater for the rectangular shape. 

23 






















































The strain and kine tic energy response between the two shap es of plates are 
shown in Figures 21, 22 and 23, for the dry, two-sides wet and one-s ide wet structures 
respectively. The rectangular plate has a higher frequency and faster damping rate. The 
rectangular plate has lo wer a mplitude of oscillations. The average energies of the 
rectangular plate are higher th an those of the square plat e. The difference in energy 
between the two plate shapes is m ore pronounc ed in the two-sides wet structure, w hich 
clearly shows the FSI effect is greates t for two-sides wet structur e and the overall FSI 
effect is less for the rectangular vice square plate since the peak energy is greater. 




Figure 21. Dry Structure Strain and Kinetic Energy Comparison for Plate Shape 
Variations with Concentrated Force and Clamped Boundary 




Figure 22. Two-sides Wet Structure Strain and Kinetic Energy Comparison for Plate 
Shape Variations with Concentrated Force and Clamped Boundary 
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Figure 23. One-side Wet Strueture Strain and Kinetie Energy Comparison for Plate 
Shape Variations with Concentrated Force and Clamped Boundary 


The comparison of strain between the two shaped plates is shown in Appendix E. 
The shear strains for the rectangular plate are all less than for the e quivalent area square 
plate. The average norm al strain at center lo cation of all three structures is a little more 
for the rectangular vice the squa re plate, while the norm al strains at the side and quarter 
locations are very sim ilar. Overall, the re duction in energy due to ESI effects of the 
rectangular plate shape is less than the square plate. 


H, COMPOSITE DENSITY 

Next the influence of com posite m aterial density on dynam ic response is 
examined. Since the response of displacement is similar to strain energy, only the strain 
energy will be used here and displacem ent plots are in Appendix F. Figures 24, 25 and 
26 show the strain and ki netic energy for the dry, two- sides wet and one-side wet 
structures respectively, w ith each u sing a composite plate of density 2020 kg/m and 
elastic m odulus of 1.7e Pa, with concentrated force and clam ped boundary for 
normalization. For the dry structu re it is clearly visible that increasing density causes a 
decrease in frequency, however, due to FSI this feature is not as pronounced in the w et 
structures. The wet structures show only s light difference in frequency and the peak 
strain energ y occurs in lowest density with on ly m inimal decrease in peak energy as 
density increases. Thekineti c energy shows a faster rate of da mping with increasing 
density for the wet structures. 
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Figure 24. Dry Structure Strain and Kinetic Energy Comparison for Density 
Variations with Concentrated Force and Clamped Boundary 




Figure 25 . Two-sides Wet Structure Strain and Kinetic Energy Comparison for 
Density Variations with Concentrated Eorce and Clamped Boundary 
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Figure 26. One-side Wet Strueture Strain and Kinetie Energy Comparison for 
Density Variations with Concentrated Force and Clamped Boundary 


To highlight the specific ESI effects, the m ethod of normalization was altered 
such that each wet structure is norm alized to its respec tive dry stru cture and the th ree 
different densities are plotted together to show its effect on response. W ith this alternate 
normalization, the strain and ki netic energy is shown in Fi gures 27 and 28 for the two- 
sides wet and one-side wet struct ures, respectively. In this representation, it is clear ESI 
gives a reduction in peak energy, is m ore signif icant in two-sid es wet structu re and 
drastically reduces the high frequency oscillat ion from the dry structure. Also, the two- 
sides wet structure kinetic energy shows faster response, as density increases, from initial 
load application to peak energy value and subsequent decay toward steady state value. 




Figure 27. Two-sides Wet Structure Strain and Kinetic Energy Comparison for 
Density Variations with Concentrated Force and Clamped Boundary (Alternate 

Normalization) 
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Figure 28. One-side Wet Strueture Strain and Kinetic Energy Comparison for 
Density Variations with Concentrated Force and Clamped Boundary (Alternate 

Normalization) 

The variations in strain of the three structures are shown in Appendix F. For the 
dry and one-side wet structure, a decreas e in frequency as d ensity increases is evid ent, 
but is not identif table f or the two- sides wets tructure. N one of the structures exhibit 
significant variation in magnitude of strain for the different density values. As shown i n 
Figure 29, using the previously discussed alternate normalization, the strains for the two- 
sides wet structure have sim ilar relative magnitude, for the three loca tions (center, side, 
quarter), with no discemable shift in frequency due to density variations, while Figure 30 
shows the one-side wet structure having a deer ease in f requency from density inc reases 
with similar relative strain magnitude. 
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Figure 29. Normal and Shear Strains for Comparison of Differ Density for Two-sides 
Wet Strueture with Concentrated Force and Clamped Boundary (Alternate 

Normalization) 
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Figure 30. Normal and Shear Strains for Comparison of Differ Density for One-side 
Wet Structure with Concentrated Force and Clamped Boundary (Alternate 

Normalization) 
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I. COMPOSITE 


MODULUS 


The influence of com posite m aterial elastic modulus on dynam ic response is 
examined next. The d isplacement response fo r different elastic m odulus is shown in 
Appendix G. Figures 31, 32 and 33 show the strain and kinetic energy for the dry, two- 
sides wet and one-side we t structures respectively, with each using a com posite plate of 
density 2020 kg/m and elastic m odulus of 1.7e Pa, with concentrated fore e and 
clamped boundary for norm alization. For the thr ee structures, an increase in frequency 
and decrease in am plitude is clearly visi ble f or increasing elastic m odulus. As the 
composite elastic modulus increases, the stru cture becomes stiffer and as the strain and 
kinetic energy plots show, the average energy decreases with increasing modulus. Also, 
the amplitude of oscillation decreases with increasing modulus. The wet structures show 
a similar rate of damping with increasing modulus. 




Figure 31. Dry Structure Strain and Kinetic Energy Comparison for Elastic Modulus 
Variations with Concentrated Force and Clamped Boundary 
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Figure 32. Two-sides Wet Structure Strain and Kinetic Energy Comparison for 
Elastic Modulus Variations with Concentrated Eorce and Clamped Boundary 



2.5 



Time (sec.) 


Eigure 33. One-side Wet Structure Strain and Kinetic Energy Comparison for Elastic 
Modulus Variations with Concentrated Eorce and Clamped Boundary 


To highlight the specific ESI effects, the m ethod of normalization was altered as 
discussed in the previous section. With this alternate normalization, the strain and kinetic 
energy is shown in Eigures 34 and 35 for the two-sides wet and one-side wet structures, 
respectively. In this representation, it is cl ear E SI gives a reduction in peak energy, is 
more significant in two-sides wet s tructure and drastically reduces the high frequency 
oscillation. The average strain energy is co mparable at each of the different elastic 
modulus values, with only the frequency and amplitude varying, while the kinetic energy 
tends to decrease with increasing modulus. 
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Figure 34. Two-sides Wet Strueture Strain and Kinetie Energy Comparison for 
Elastic Modulus Variations with Concentrated Eorce and Clamped Boundary 

(Alternate Normalization) 




Eigure 35. One-side Wet Structure Strain and Kinetic Energy Comparison for Elastic 
Modulus Variations with Concentrated Eorce and Clamped Boundary (Alternate 

Normalization) 

The strain variations of the three structures are shown in Appendix G. Eor these 
structures, the decrease in amplitude of oscillation as elastic modulus increases is evident 
while the in crease in f requency les s noticeable in the strain plots. The m ost notable 
feature is the large re duction in strain as the m odulus increases; this is due to the 
increased stiffness. Us ing the alternate norm alization, so that spec ific influences of 
elastic modulus and ES I can be highlighted fo r strain, the strain responses of the two- 
sides wet and one-side wet structures are shown in Eigures 36 and 37 , respectively. Eor 
the two-sides wet structure, th e relative magnitude of strain is consistent for each of the 
three moduli across the three different locations on the plate, with minor indication of the 
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increase in frequency with increasing m odulus. The frequency shift is m uch more 
evident in the one -side wet s tructure, whil e again the rela tive stra in m agnitude is 
consistent for the various modulus values. 







Figure 36. Normal and Shear Strains for Comparison of Differ Elastic Modulus for 
Two-sides Wet Structure with Concentrated Foree and Clamped Boundary 

(Alternate Normalization) 
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Figure 37. Normal and Shear Strains for Comparison of Differ Elastic Modulus for 
One-side Wet Structure with Concentrated Force and Clamped Boundary 

(Alternate Normalization) 
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J. IMPACT 


LOADING 


The final dynam ic behavior exam ined is the impact response of com posite plate 
from a steel projectile. Three velocities are examined for the projectiles: 1 m/s, 5 m/s and 
10 m/s. In addition, the response due to tw o impact face shapes are com pared, each 
having the same surface area for i mpact and equal mass. A cylindrical shaped impacto r 
has a circular shape area of impact and a rectangular shaped impactor has a square shape 
area of impact. Each projectile contacts the composite plate at the center. 

1, Shape of Impactor 

To investigate any dependence on shape of the impact object, a fixed velocity of 
10 m/s is used to com pare the difference in response between the tw o shapes of im pact 
projectiles. Figures 38, 39 and 40 compare the displacement, strain energy and kinetic 
energy response of the three structu res to the circular and square shape impactor. The 
figures show ESI gives a significant reduction in amplitude and frequency and the square 
impact face has less am plitude for dry and one-side wet s tructure than the circula r face 
impactor. The two-sides wet structure show s similar initial response between the two 
shapes of impact and then the square face impactor respo nse stead ies out with higher 
amplitude o f oscillation . The average di splacement is co mparable b etween th e two 
impact shapes, while the average energy is less for the square face impact for each of the 
three structures. 




Circular Impact Area 


Square Impact Area 


Figure 38. Displacement Comparison of Three Structures Due to Different Impactor 

Shapes at 10 m/s 
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Circular Impact Area Square Impact Area 

Figure 39. Strain Energy Comparison of Three Struetures Due to Different Impaetor 

Shapes at 10 m/s 
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Circular Impact Area Square Impact Area 

Figure 40. Kinetic Energy Comparison of Three Structures Due to Different Impaetor 

Shapes at 10 m/s 


The norm al strain at each lo cation is compared in Figure 41 for the two 
shapes of impaetor. The ESI decreas es si ightly th e str ain amplitude with m inor 
decrease in frequency. The two-sides wet st ructure has a larger peak strain for the 
center and side locations for th e square faced im pactor over the circular faced 
impaetor, while the cylindrical impaetor ha s peak strain in dry and one-side wet 
structure. The shear strain is slightly higher for square face im pactor at center 
location, while the shear strains are com parable between im pactor at the side and 
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quarter locations. The FSI is m ore pronounced for the square impactor one-side 
wet, and more for the cylindrical impactor two-sides wet case. 




2.5 3 3.5 

Time (sec.) 






Circular Impact Area 


Square Impact Area 


Figure 41. Normal Strain Comparison of Three Structures Due to Different Impactor 

Shapes at 10 m/s 
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Circular Impact Area 


Square Impact Area 


Figure 42. Shear Strain Comparison of Three Structures Due to Different Impactor 

Shapes at 10 m/s 
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Additional figures comparing the response between the square and circular faced 
impactor are contained in Appendix H. These additional figures highlight the FSI 
differences for the two different shapes of impactor for the two-sides wet and one-side 
wet structures by norm alizing each to their respective dr y structure response. To 
summarize, FSI slightly decreases the strain amplitude and frequency of oscillation with 
the square impactor having a slight increase in frequency and am plitude of oscillation 
over the circular impactor for both two-sides and one-side wet. The square impactor also 
has slightly less av erage energy (strain and kinetic). The stra ins are nearly the s ame 
between the cylinder and square im pactor with the square having higher peak strain at 
center position and comparable for the side and quarter positions. The average strains are 
roughly the same except for the center pos ition which is higher due to higher peak strain 
initially. 

2. Velocity of Impact 

The effect of impact velocity is straight forward; increasing impact velocity gives 
increased magnitude of plate displacement, strain and kinetic energies. When combining 
the varying im pact velocities with different s haped im pactors, there are som e slight 
differences in response. The shift in initial response when comparing the three velocities 
is due to th e time difference required for the impactor to trave rse the distance to the 
composite plate and should not be misinterpreted as a frequency shift. 

The response of each o f three structures to different initial im pact velocities is 
shown in A ppendix I for both circular and square faced impactors. Increas ing impact 
velocity s imply increas es the respon se. Ge nerally, the square faced impactor has less 
amplitude of oscillation and average values for the dry and one-side wet structure and the 
two-sides wet structure amplitude of oscillation and average value is s imilar for the two 
different impactors. The normal strains are comparable with only very slight decrease in 
amplitude o f oscillation for the square impactor for each o f the th ree structures. The 
shear strain is also similar am ong the three structures and two im pactors for the three 
impact velocities. 
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To focus on the FSI effects, each o f the impact velocities for the two-s ides and 
one-side wet stru ctures are norm alized to a respective dry structure. T hese normalized 
responses are shown in Figures 43, 44 and 45 fo r displacement, strain energy and kinetic 
energy respectively. Using this normalization, it is clear FSI causes significant decreases 
in frequency and a mplitude range of res ponse. The two-sides wet structure shows 
decreased peak values while the on e-side wet has slightly increased peak values over the 
strictly dry structure du e to effects of the wa ter layer on one side of plate. The square 
faced impactor has higher relative amplitude of oscillation for displacement, but slightly 
lower energies than those of circular faced impactor. 




Time (sec.) 




Circular Square 

Figure 43. Comparison of Displacement Response Due to Impact Velocity Effects for 

Circular and Square Faced Impactor 
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Circular Square 

Figure 44. Comparison of Strain Energy Response Due to Impact Velocity Effects for 

Circular and Square Faced Impactor 


The normal and shear strains using the nor malization to highlight the F SI effects 
are shown in Figures 46 and 47, respectively. These figures show the relative magnitudes 
of strain f or all three velocitie s ar e sim ilar with the ex ception of the two-sides wet 
structure with square face impactor has slightly increased peak strain. 
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Cireular Square 

Figure 45. Comparison of Kinetie Energy Response Due to Impaet Veloeity Effects 

for Circular and Square Eaced Impactor 
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Circular Square 

Figure 46. Comparison of Normal Strain Due to Impact Velocity Effects for Circular 

and Square Faced Impactor 
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Figure 47. Comparison of Shear Strain Due to Impact Velocity Effects for Circular 

and Square Faced Impactor 
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V. NUMERICAL MODEL COMPARISON TO EXPERIMENT 


The principal focus of this work is nu merical study of various param eters which 
affect the dynam ic be havior of com posite plates. A separate study conducted 
experimentally exam ines the behavior of dry and wet plates s ubjected to impact. 
Numerical and experim ental studies each have their respective advantages and 
disadvantages and are used to com plement each other. In particular, the v arious 
parametric studies conducted in this wo rk were only possible utilizing numerical 
modeling. Experimental testing is limited to measuring forces and strains through gages. 
Preliminary comparison of experim ental and nu merical work is noted here to determine 
methods for improvement to follow on research. 

A, EXPERIMENTAL SETUP FOR IMPACT LOADING 

Preliminary experimental behavior study of thin composite plates is conducted on 
12 inch square and 1/16 inch thick plate clamped in the frame of an impact testing device. 
The device uses a weighted sled system to strike a cylindrical impactor. Multiple impacts 
are prevented using a large sp ring opposing the cylindrical impactor such that only one 
impact event takes place. There is a for ce measuring gage m ounted on the end of the 
cylindrical impactor to m easure the force durin g contact with the com posite plate. The 
schematic of the experimental device setup is shown in Figure 48. The underside of the 
composite plate is instrumented with strain gages, bonded to the pi ate with epoxy, in the 
layout shown in Figure 49. Th e strain gages measure approxi mately 1 cm square. Gage 
2 is in the center location for comparison to the numerical model a nd is directly below 
impact site. Gages 1 and 4 are repr esentative of a side location sim ilar to the num erical 
model and gages 3 an d 5 are sim ilar to the quarter location. The data acqu isition 
software measures the transient force and strain data at 1000 Hz sampling rate. 

B, NUMERI CAL MODEL 

The 1/16**^ inch thick, 12 inch by 12 inch com posite plate is m odeled using shell 

elements with a m esh seed of 60 nodes per side . The m esh size is chosen to adequately 

approximate the im pact force gage area and strain gage size reasonab ly. The i mpact 
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force m easured experim entally is converted to equivalent pressure and applied to 
elements approximating a cylindrical impactor striking the plate. In the numerical model, 
stress is co mputed and strains are calcu lated using standard stress-strain transform ation 
equations. The strain over the area of the numbered experim ental strain gages is 
calculated by averaging the elem ents which ap proximate the size of the strain gage to 
compare with the experimentally measured strains. 


Impuct Icstinj; machine 


Holder for 

crane lifting O 


Sliding — 
Impact weight 


Coil Spring 


Supporter for 
hanging the frame- 


Guidcr / 
for Rad ^— 


Impact — 
transmitting rod 


Caiiipositc plate' 


n 


P ' v| 


48 in. 


5 in. 


12 in. 


12 in. 


12 in. 


5 in. 



(a) Schematic 


(b) As-built 


Figure 48. Impact Device Experimental Setup 
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(a) Strain Gage Layout (b) Actual Composite Plate 

Figure 49. Experiment Strain Gage Layout on Underside of Composite Plate 
(Dimensions in parenthesis are given in inches) 


C. COMPARISON OF EXPERIMENTAL AND NUMERICAL RESULTS 

One of the challenges experim ental work presents when m ensuring strains is 
getting a good bond between the strain gages an d the composite plate. Another is having 
the strain g ages aligne d perf ectly with the direction of fibers in the com posite and 
minimizing the area of the gage covering th e matrix which f orms the com posite. The 
experiment is conducted in a one-side wet sc enario using an anechoic tank to m inimize 
water disturbance effects, with the side oppos ite of i mpact on the composite plate kept 
dry through a plexi-glass box bound to the underside of the com posite plate. The 
experiment is also run in a com pletely dry condition. Both dry and wet cases use the 
same impact force by dropping the weighted sled from full height, giving the steel impact 
rod roughly a 5 m /s initial ve locity. The experim ents were first conducted in the wet 
condition and then dry. Following the wet experim ents it was identified that the strain 
gage labeled Gage 1 had broken free from the com posite plate and hence was not 
available for the dry experiment. 

The comparison of normal strain for gage 1 location, between the experim ent and 
simplified Finite Element model is shown in Figure 50. As shown, there is not a good 
comparison between the experim ent and m odel results for the one-side wet condition 
(note that the strain gage fell off prior to dry experiment and there is no experimental data 
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to compare). Because the strain gage fell off after the wet experiment, the data shown is 
possibly erroneous due to the strain gage disbonding and any com parison at the gage 1 
location is suspect. 






Dry Condition 


One-side Wet Condition 


Figure 50. Comparison of Normal Strain at Gage 1 Location Between Experiment 

and FEM in Dry and One-side Wet Condition 


The comparison of normal strain for gage 2 location, between the experim ent and 
simplified FEM m odel is shown in Figure 51. As shown, the com parisons between the 
experiment and model results are quite good for both the one-side wet and dry condition. 
This good agreem ent between the experim ent and numerical model is evidence of the 
feasibility to aceurately pred ict co mposite plate response us ing finite elem ent models. 
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This gives more flexibility for researehers as many more parameters ean be varied with a 
numerieal model. The faet that th e x-axis model strain is higher than the experim ent and 
the y-axis is lower is an indication there m ay be som e misalignment of t he strain gage 
with the fiber direction. If this is the case, som e improvement can be obtained through 
use of a Mohr Circle transformation. 






Dry Condition 


One-side Wet Condition 


Figure 51. Comparison of Normal Strain at Gage 2 Location Between Experiment 

and FEM in Dry and One-side Wet Condition 


The comparison of normal strain for the gage 3 location, between the experim ent 
and simplified FEM model is shown in Figure 52. As shown, the com parisons between 
the experiment and model results are quite good for the dry condition but not the one-side 
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wet condition. Use of Mohr Circle transf ormation m ay i mprove the dry com parison. 
What is encouraging is the trend between e xperiment and model tr acks. Unfortunately, 
there is not a good explanation of why the dr y condition is in such good agreem ent but 
the one-side wet condition is not. 






Dry Condition 


One-side Wet Condition 


Figure 52. Comparison of Normal Strain at Gage 3 Location Between Experiment 

and FEM in Dry and One-side Wet Condition 


The comparison of normal strain for gage 4 location, between the experim ent and 
simplified FEM m odel is shown in Figure 53. As shown, there is good agreem ent 
between the experiment and model results for the dry cond ition and the trends agree for 
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the one-side wet condition, although the m agnitudes are off. Again a Mohr Circle 
transformation could improve the dry and wet comparison. 






Dry Condition 


One-side Wet Condition 


Figure 53. Comparison of Normal Strain at Gage 4 Location Between Experiment 

and FEM in Dry and One-side Wet Condition 


The comparison of normal strain for gage 5 location, between the experim ent and 
simplified FEM m odel is shown in Figure 54. As shown, the trends between the 
experiment and m odel are sim ilar, but the magnitudes are not, and application of Mohr 
Circle will not im prove the va lues as both the x-axis and y- axis normal strains are over 
predicted in the numerical model. 
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Dry Condition 


One-side Wet Condition 


Figure 54. Comparison of Normal Strain at Gage 5 Location Between Experiment 

and FEM in Dry and One-side Wet Condition 


In summary, the strain in the vicinity of impact, at gage 2 location, compares very 
well between model and experiment. Moving away from impact location either to side or 
quarter lo cation r esults in le ss agr cement be tween the experim ent and m odd. This 
indicates that proper strain gage alignm ent with fiber direction and good bonding over 
fiber vice matr ix is im portant. Other things to consider in f uture work are altering the 
element size in the m odd and adjusting the quantity of elements used in aver aging to 
determine the strain at a gage location for comparison to the experimental data. 
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VI. CONCLUDING REMARKS AND RECOMMENDATIONS 


Thin com posite p late structures were exam ined under various con ditions to 
investigate the effect of FSI on dynam ic be haviors. Overall, water influenced 
significantly both kinetic and st rain energies of the com posite structures by greatly 
reducing their m agnitudes and frequencies. The FSI greatly suppressed the oscillatory 
nature of dynamic responses of the structures . Whether a structure is wet on one-side or 
two-sides, the FSI effect was very clear even though the two-sides wet structures showed 
a greater FSI effect. 

The boundary condition, either clamped or simple, has sim ilar behaviors and is 
thereby not a significant contributor for FSI. The size and sh ape of the composite plates 
was shown to have m inor differences in FSI. The m ethod of loading the plate, either 
concentrated force, uniform pressure or impact, showed some difference on the degree of 
FSI. Interestingly, the shape of t he im pacting object (contact shape) gave different 

degrees of FSI for equivalent im pact velocities. The larg est variation of FSI was due to 
differences in m aterial properties such as de nsity and elastic m odulus. As a result, it is 
critical to understand and incorporate the FSI effects when designing reliable composite 
structures employed in an underwater environment. 

Future work should exam ine the dynam ic behavior of composites which include 
moisture ab sorption ef fects. Additiona lly, various types of co mposites should be 
compared for determ ination of th e best response behavior properties and m inimal 
moisture absorption. Finally, both num erical and experim ental work should be 
conducted to m onitor com posite behavior in failure. The failure modes should be 
investigated as to wheth er they are m atrix or fiber failure, delam ination or a m ixture of 
failure modes. 
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Norm. Displacement, CntrTop Plate, 1Wet Norm. Displacement, CntrTop Plate, 2Wet Norm. Displacement, CntrTop Plate, Dry 


APPENDIX A: ADDITIONAL FIGURES FOR CLAMPED AND 
SIMPLE BOUNDARY WITH CONCENTRATED FORCE LOAD 








Figure 55. Displacement and Strain Energy Comparison of Clamped and Simple 

Boundary with Concentrated Force Load 
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Figure 56. Normal and Shear Strain Comparison at Center Position for Clamped 
versus Simple Boundary with Concentrated Force Load 
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Norm.e ,Side Bottom Plate, 1 Wet Norm.e ,Side Bottom Plate, 2Wet Norm.e , Side Bottom Plate, Dry 


Time (sec.) 



Time (sec.) 






Figure 57. Normal and Shear Strain Comparison at Side Position for Clamped versus 

Simple Boundary with Concentrated Force Load 
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APPENDIX B: ADDITIONAL FIGURES FOR FORCE AND 
PRESSURE LOAD COMPARISON WITH CLAMPED BOUNDARY 





Figure 58. Comparison of Dry Structure Response for Displacement, Strain and 
Kinetic Energies Between Force and Pressure Loading with Clamped Boundary 
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Norm. , Qrtr Bottom Plate, Dry Norm, , Side Bottom Plate, Dry Norm, Cntr Bottom Plate, Dry 







Figure 59. Normal and Shear Strains for Comparison of Dry Structure with Clamped 

Boundary between Force and Pressure Loading 
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Figure 60. Normal and Shear Strains for Comparison of Two-sides Wet Structure 
with Clamped Boundary between Force and Pressure Loading 
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Norm.e , Qrtr Bottom Plate, 1 Wet Norm. 8 , Side Bottom Plate, 1 Wet Norm.e , Cntr Bottom Plate, 1 Wet 








Figure 61. 


Normal and Shear Strains for Comparison of One-side Wet Strueture with 
Clamped Boundary between Foree and Pressure Loading 
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APPENDIX C: ADDITIONAL FIGURES FOR FORCE AND 
PRESSURE LOAD COMPARISON WITH SIMPLE BOUNDARY 





Figure 62. Comparison of Dry Structure Response for Displacement, Strain and 
Kinetic Energies between Force and Pressure Loading with Simple Boundary 
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Figure 63. 


Normal and Shear Strains for Comparison of Dry Structure with Simple 
Boundary between Force and Pressure Loading 
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Figure 64. Normal and Shear Strains for Comparison of Two-sides Wet Structure 
with Simple Boundary between Force and Pressure Loading 
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Figure 65. Normal and Shear Strains for Comparison of One-side Wet Structure with 
Simple Boundary between Force and Pressure Loading 
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APPENDIX D: ADDITIONAL FIGURES FOR PLATE SIZE 
EFFECTS WITH CONCENTRATED FORCE LOAD AND 
CLAMPED BOUNDARY 





Figure 66. Comparison of Displacement Response for Three Structures Due to Size 
Variation Effects with Concentrated Force and Clamped Boundary 
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Norm. 6^, Qrtr Bottom Plate, Dry Norm. , Side Bottom Plate, Dry Norm.e , Cntr Bottom Plate, Dry 








Figure 67. Normal and Shear Strains for Comparison of Differ Plate Sizes for Dry 
Structure with Concentrated Force and Clamped Boundary 
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Norm. , Qrtr Bottom Plate, 2Wet Norm, e^ , Side Bottom Plate, 2Wet Norm.E , Cntr Bottom Plate, 2Wet 








Figure 68. Normal and Shear Strains for Comparison of Differ Plate Sizes for Two- 
sides Wet Strueture with Concentrated Force and Clamped Boundary 
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Norm.Ej^.Qrtr Bottom Plate, 1 Wet Norm.,Side Bottom Plate, 1 Wet Norm.e , Cntr Bottom Plate, 1 Wet 








Figure 69. Normal and Shear Strains for Comparison of Differ Plate Sizes for One- 
side Wet Strueture with Concentrated Force and Clamped Boundary 
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APPENDIX E: ADDITIONAL FIGURES FOR PLATE SHAPE 
EFFECTS WITH CONCENTRATED FORCE LOAD AND 
CLAMPED BOUNDARY 




Time (sec.) 



Figure 70. Comparison of Displacement Response for Three Structures Due to Shape 
Effects with Concentrated Force and Clamped Boundary 
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Figure 71. Normal and Shear Strains for Comparison of Differ Plate Shapes for Dry 
Strueture with Concentrated Force and Clamped Boundary 
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Norm. E^.Qrtr Bottom Plate, 2Wet Norm., Side Bottom Plate, 2Wet Norm.E , Cntr Bottom Plate, 2Wet 







Figure 72. 


Normal and Shear Strains for Comparison of Differ Plate Shapes for Two- 
sides Wet Structure with Concentrated Force and Clamped Boundary 
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Norm.E^, Qrtr Bottom Plate, 1 Wet Norm. e^, Side Bottom Plate, t Wet Norm.E , Cntr Bottom Plate, tWet 








Figure 73. Normal and Shear Strains for Comparison of Differ Plate Shapes for One- 
side Wet Structure with Concentrated Force and Clamped Boundary 
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APPENDIX F: ADDITIONAL FIGURES FOR COMPOSITE 
DENSITY EFFECTS WITH CONCENTRATED FORCE LOAD AND 

CLAMPED BOUNDARY 


The following use composite density of 2020 kg/m ^ and modulus 1.7e GPa for 
normalization in each of the three structures. 


a) Dry Structure 





Figure 74. Comparison of Displacement Response for Three Structures Due to 
Density Effects with Concentrated Force and Clamped Boundary 
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Figure 75. Normal and Shear Strains for Comparison of Differ Density for Dry 
Structure with Concentrated Force and Clamped Boundary 
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Norm. ,Qrtr Bot. Plate, 2Wet Norm, e^ .Side Bot. Plate, 2Wet Norm.E ,Cntr Bot. Plate, 2Wet 
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Figure 76. Normal and Shear Strains for Comparison of Differ Density for Two-sides 
Wet Strueture with Concentrated Force and Clamped Boundary 
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Norm. 6^ .QrtrBot. Plate, 1 Wet Norm, .Side Bet. Plate, t Wet Norm.s ,Cntr Bot. Plate, t Wet 
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Figure 77. Normal and Shear Strains for Comparison of Differ Density for One-side 
Wet Strueture with Concentrated Force and Clamped Boundary 
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APPENDIX G: ADDITIONAL FIGURES FOR COMPOSITE 
ELASTIC MODULUS EFFECTS WITH CONCENTRATED FORCE 
LOAD AND CLAMPED BOUNDARY 


The following use composite density of 2020 kg/m ^ and modulus 1.7e GPa for 
normalization in each of the three structures. 


a) Dry Structure 



b) Two-sides wet Structure 



c) One-side Wet Structure 



Figure 78. Comparison of Displacement Response for Three Structures Due to 
Elastic Modulus Effects with Concentrated Force and Clamped Boundary 
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Figure 79. Normal and Shear Strains for Comparison of Differ Elastie Modulus for 
Dry Structure with Concentrated Force and Clamped Boundary 
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Figure 80. Normal and Shear Strains for Comparison of Differ Elastie Modulus for 
Two-sides Wet Structure with Concentrated Force and Clamped Boundary 
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Norm. e^.QrtrBot. Plate, 1 Wet Norm, .Side Bot. Plate, t Wet Norm.e .CntrBot. Plate, t Wet 








Figure 81. Normal and Shear Strains for Comparison of Differ Elastie Modulus for 
One-side Wet Structure with Concentrated Force and Clamped Boundary 
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APPENDIX H: ADDITIONAL FIGURES FOR IMPACTOR SHAPRE 
EFFECTS WITH CLAMPED BOUNDARY 


The following com pare circular face to s quare face im pactor, with equal im pact 
area and equal m ass, for two-sides and one-sid e wet structures norm alized to respectiv e 
dry structure. 



Two-sides Wet One-side Wet 


Figure 82. Comparison of Displacement Response for Two-sides and One-side Wet 

Structures Due to Impactor Shape Effects 
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Figure 83. Comparison of Strain Energy Response for Two-sides and One-side Wet 

Structures Due to Impactor Shape Effects 
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Figure 84. Comparison of Kinetic Energy Response for Two-sides and One-side Wet 

Structures Due to Impactor Shape Effects 
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Eigure 85. Normal and Shear Strain Comparison of Different Impactor Shape for 
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Two-sides and One-side Wet Structure at Center Location 
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Figure 86. Normal and Shear Strain Comparison of Different Impactor Shape for 
Two-sides and One-side Wet Structure at Side Location 
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Figure 87. Normal and Shear Strain Comparison of Different Impactor Shape for 
Two-sides and One-side Wet Structure at Quarter Location 
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APPENDIX I: ADDITIONAL FIGURES FOR IMPACT VELOCITY 

AND SHAPE EFFECTS 


The following com pare circular face to s quare face im pactor, with equal im pact 
area and equal mass, for different velocities of the three structures, normalized to 1 m/s. 


Dry 




Time (sec.) 






Circular Square 

Comparison of Displacement Response for Three Structures Due to 
Impact Velocity Effects for Circular and Square Faced Impactor 
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Figure 88. 
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Figure 89. Comparison of Strain Energy Response for Three Struetures Due to 
Impact Velocity Effects for Circular and Square Eaced Impactor 
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Figure 90. Comparison of Kinetie Energy Response for Three Structures Due to 
Impact Velocity Effects for Circular and Square Eaced Impactor 
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Figure 91. Comparison of Normal Strain at Center Location for Three Structures Due 
to Impact Velocity Effects for Circular and Square Faced Impactor 


92 











































































































































































8 


Dry 



.4 1_^^^_I_I_^_I_I_^_I 

0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 O.C 

Time (sec.) 


0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 O.i 

Time (sec.) 


Two 


sides 

wet 




One 


side 

wet 




Cireular 


Square 


Figure 92. Comparison of Shear Strain at Center Loeation for Three Structures Due 
to Impact Velocity Effects for Circular and Square Faced Impactor 
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